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Topological  insulators  exhibit  metallic  surface  states  populated  by  massless  Dirac  fermions  with  spin-momentum  locking, 
where  the  carrier  spin  lies  in-plane,  locked  at  right  angles  to  the  carrier  momentum.  Here,  we  show  that  a  charge  current 
produces  a  net  spin  polarization  via  spin-momentum  locking  in  Bi2Se3  films,  and  this  polarization  is  directly  manifested  as 
a  voltage  on  a  ferromagnetic  contact.  This  voltage  is  proportional  to  the  projection  of  the  spin  polarization  onto  the 
contact  magnetization,  is  determined  by  the  direction  and  magnitude  of  the  charge  current,  scales  inversely  with  Bi2Se3 
film  thickness,  and  its  sign  is  that  expected  from  spin-momentum  locking  rather  than  Rashba  effects.  Similar  data  are 
obtained  for  two  different  ferromagnetic  contacts,  demonstrating  that  these  behaviours  are  independent  of  the  details  of 
the  ferromagnetic  contact.  These  results  demonstrate  direct  electrical  access  to  the  topological  insulators'  surface-state 
spin  system  and  enable  utilization  of  its  remarkable  properties  for  future  technological  applications. 


Topological  insulators  (TIs)  form  a  new  quantum  phase  of 
matter  distinct  from  the  classic  dichotomy  of  metals  and  semi¬ 
conductors1-4.  Whereas  the  bulk  states  form  a  bandgap,  the 
surface  states  form  a  Dirac  cone  similar  to  graphene  (Fig.  la)  and 
are  topologically  protected  against  disorder  scattering.  In  marked 
contrast  with  the  spin-degenerate  bands  of  graphene,  TI  surface 
states  are  spin  polarized.  Topological  insulators  are  expected  to 
produce  new  functionalities5  and  enable  insights  into  complex 
phenomena  in  many  scientific  arenas  including  spintronics6-8, 
quantum  information  technology9,  highly  correlated  electron 
systems10,  magnetic  monopoles11  and  quantized  magnetoelectric 
coupling12.  Examples  of  TI  materials  include  Bij^Sb^  (ref.  4), 
Bi2Se3,  Bi2Te3  and  Sb2Te3  (refs  13-15). 

One  of  the  most  striking  properties  is  spin-momentum  locking; 
the  spin  of  the  TI  Dirac  surface  state  lies  in-plane,  and  is  locked  at 
right  angles  to  the  carrier  momentum.  An  unpolarized  charge 
current  should  thus  create  a  net  spin  polarization,  the  amplitude 
and  orientation  of  which  are  controlled  by  the  charge  current. 
This  remarkable  property  has  been  anticipated  by  theory16-18,  but 
never  accessed  in  a  simple  transport  structure.  A  top  view  of  a 
cross-section  of  the  surface-state  bands  taken  above  the  Dirac 
point  is  approximately  circular  (with  slight  hexagonal  warping19), 
with  the  spin  tangential  at  all  points,  as  shown  by  the  solid  line 
in  Fig.  lb.  A  net  momentum  along  the  kx  direction  (represented 
by  a  displacement  of  this  circular  cross-section  along  kx,  indicated 
by  the  dashed  circle)  produced  by  an  electric  field  results  in  an 
electron  current  7e  along  x.  Owing  to  spin-momentum  locking, 
this  simultaneously  induces  a  spin  current  with  a  spin  polarization 
oriented  along  —y.  Thus,  an  unpolarized  surface-state  charge 
current  creates  a  net  spin  polarization,  with  amplitude  and  orien¬ 
tation  determined  by  the  amplitude  and  direction  of  the  charge 
current  (Fig.  lc).  Spin-momentum  locking  has  been  probed 
by  spin-resolved  photoemission14'20'21  and  polarized  optical 
spectroscopic  techniques22,  but  never  accessed  in  a  simple 
transport  geometry. 


Figure  1  |  Schematic  of  TI  surface  bands  and  experimental  concept. 

a,  Dirac  cone  of  the  TI  surface  states  (blue),  with  the  spin  at  right  angles  to 
the  momentum  at  each  point.  The  bulk  conduction  and  valence  bands  are 
shown  in  grey,  b,  Top  view  of  a  slice  in  the  kx—ky  plane  of  the  TI  surface 
states.  An  applied  voltage  produces  a  net  momentum  along  kx  and  spin- 
momentum  locking  gives  rise  to  a  net  spin  polarization  oriented  in-plane  and 
at  right  angles  to  the  current,  c,  Concept  drawing  of  the  transport 
experiment.  The  voltage  measured  at  the  ferromagnetic  detector  is 
proportional  to  the  projection  of  the  current-induced  TI  spin  polarization 
onto  the  contact  magnetization. 
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Figure  2  |  Schematic  of  contacts  and  transport  devices,  a,  Cross-section  of  a  typical  ferromagnetic  contact,  b,  Top  view  of  contact  layout  for  devices  with  a 
single  row  of  collinear  ferromagnetic  detector  contacts  (solid  red)  and  coordinate  axes.  The  ferromagnetic  contacts  have  dimensions  of  80  x  80  pm 
separated  by  45  pm,  edge  to  edge.  The  representative  cross-section  does  not  show  all  of  the  contacts,  for  clarity,  c,  As  in  b,  but  for  devices  with  two  parallel 
rows  of  collinear  detector  contacts:  top  row  is  ferromagnetic  (solid  red);  bottom  row  is  a  non-magnetic  reference.  The  ferromagnetic  contacts  have 
dimensions  of  either  80  x  80  pm  or  60  x  60  pm,  separated  by  45  pm  and  40  pm  edge  to  edge,  respectively.  An  adjacent  parallel  row  of  non-magnetic 
reference  contacts  (Au/Ti)  of  the  same  dimensions  is  separated  from  the  ferromagnetic  contacts  by  15  pm,  edge  to  edge.  The  representative  cross-section 
does  not  show  all  of  the  contacts  for  clarity. 


Here,  we  show  direct  electrical  detection  of  the  bias  current- 
induced  spin  polarization  as  a  voltage  measured  on  a  ferromagnetic 
metal  tunnel  barrier  surface  contact.  The  magnetization  of  the 
contact  determines  the  spin  detection  axis,  and  the  voltage  measured 
at  this  contact  is  proportional  to  the  projection  of  the  TI  spin  polar¬ 
ization  onto  this  axis.  Although  high  bulk  carrier  densities  provide  a 
parallel  bulk  conduction  channel,  we  observe  dependencies  on  exper¬ 
imental  parameters  that  are  most  straightforwardly  explained  by  spin- 
momentum  locking  in  surface  states.  An  unpolarized  bias  current  is 
applied  between  two  non-magnetic  contacts.  When  the  charge 
current  is  orthogonal  to  the  magnetization  of  the  ferromagnetic 
detector  contact,  the  TI  spin  is  parallel  (or  antiparallel)  to  the  magne¬ 
tization,  as  shown  in  Fig.  lc,  and  a  spin- related  signal  is  detected  at 
the  ferromagnetic  contact  proportional  to  the  magnitude  of  the 
charge  current.  When  the  direction  of  the  charge  current  is  reversed, 
the  measured  voltage  changes  sign.  As  the  contact  magnetization  is 
reversed  by  applying  a  small  in-plane  magnetic  field,  the  measured 
spin  voltage  mirrors  the  hysteresis  loop  of  the  contact.  When  the 
contact  magnetization  is  rotated  in-plane  90°  so  that  the  charge 
current  is  parallel  to  the  magnetization,  no  spin  voltage  is  detected, 
because  the  TI  spin  polarization  is  now  orthogonal  to  the  contact 
magnetization.  The  voltage  measured  scales  inversely  with  Bi2Se3 
film  thickness  and  its  sign  is  that  expected  from  spin-momentum 
locking  and  opposite  that  of  a  Rashba  effect.  Similar  data  are  obtained 
for  two  different  ferromagnetic  contact  structures,  Fe/Al203  and 
Co/MgO/graphene,  underscoring  the  fact  that  these  behaviours  are 
due  to  bias  current-induced  spin  polarization  in  the  TI  surface 
states  rather  than  the  bulk  and  are  independent  of  the  details  of  the 
contact.  These  results  demonstrate  simple  and  direct  electrical 


access  to  the  TI  Dirac  surface-state  spin  system  and  provide  clear 
evidence  for  the  spin-momentum  locking  and  bias  current-induced 
spin  polarization.  These  results  also  identify  a  pathway  for 
future  experiments  and  enable  utilization  of  TI  properties  for 
technological  applications5-8. 

Bi2Se3  is  regarded  as  a  model  TI.  It  exhibits  a  single  Dirac  cone 
and  the  largest  bandgap  (~0.3  eV)  of  known  TIs13,  and  its  topolo¬ 
gical  surface  states  have  been  shown  to  persist  to  at  least  room  temp¬ 
erature  due  to  weak  electron-phonon  coupling14,23.  Recent  work  has 
shown  that  these  topologically  protected  surface  states  are  remark¬ 
ably  insensitive  to  both  magnetic  and  non-magnetic  surface  adsor¬ 
bates24,25  and  maintain  their  linear  dispersion  even  when  forming  a 
planar  interface  with  Fe7Se8,  a  ferromagnet  with  in-plane  magneti¬ 
zation26.  Those  authors  also  found  little  solubility  of  Fe  into  the 
Bi2Se3  (ref.  26),  and  the  studies  did  not  address  spin  polarization. 
Spin-resolved  photoemission  on  ferromagnetic  ternary  alloys  with 
a  magnetic  dopant,  Bi2_xMnxSe3  (x<0.1),  showed  dramatic 
effects  on  spin  texture  near  the  Dirac  point,  but  little  at  higher  ener¬ 
gies27.  Such  robustness  is  essential  for  the  realization  of  future  device 
structures  based  on  the  interplay  of  spin  transport,  ferromagnetism 
and  spin-momentum  locking. 

Bi2Se3(lll)  films  11-45  nm  thick  were  grown  by  molecular 
beam  epitaxy  (MBE)  on  epitaxial  graphene  (EG)/SiC(0001)  and 
Al203(0001)  at  substrate  temperatures  of  275-325  °C  (see 
Methods).  The  as-grown  films  were  n-type  with  a  net  carrier 
density  ofn^lxlO19  cm  3  so  that  the  Fermi  level  was  within 
the  bulk  conduction  band. 

Two  types  of  ferromagnetic  detector  contact  were  used — 
Fe/Al203  and  Co/MgO/graphene  (see  Methods) — and  similar 
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Figure  3  |  Tl  spin  polarization  detected  as  a  voltage  with  Fe/Al203  contacts.  The  contact  layout  of  Fig.  2b  is  used  and  the  Bi2Se3  film  is  45  nm  thick. 
a,c,  Field  dependence  of  the  voltage  measured  using  contact  Fe2  as  the  detector  with  the  magnetization  collinear  with  the  induced  Tl  spin  for  bias  currents  of 
+2  mA  and  -2  mA,  respectively.  b,d,  Similar  data  for  contact  Fel  as  the  detector.  e,f,  Field  dependence  of  the  voltage  measured  using  contact  Fel  as  the 
detector  with  the  magnetization  orthogonal  to  the  induced  Tl  spin  for  bias  currents  of  +2  mA  and  —2  mA,  respectively.  Black  traces  correspond  to  positive 
magnetic  field  sweeps  (from  negative  to  positive  values)  and  red  traces  correspond  to  negative  field  sweeps.  The  relative  orientations  of  le,  induced  Tl  spin 
polarization  s  and  detector  magnetization  M  are  shown  in  the  insets.  All  data  were  obtained  at  T=  8  K. 


results  were  obtained  for  each.  Such  contacts  have  successfully 
enabled  electrical  detection  of  spin  current  and  accumulation  in 
other  material  systems28  and  their  application  to  TIs  has  recently 
been  treated  theoretically29.  The  oxide  layers  ensure  that  no  direct 
magnetic  exchange  interaction  occurs  between  the  ferromagnetic 
metal  and  Tl  surface  states.  The  samples  were  processed  into 
the  device  structures  presented  in  Fig.  2,  which  consist  of  a  co- 
linear  series  of  ferromagnetic  detector  and  non-magnetic  Au/Ti 
reference  contacts,  with  adjacent  contacts  separated  by 
40-45  |xm.  The  magnetic  characteristics  of  the  ferromagnetic  con¬ 
tacts  were  determined  from  anisotropic  magnetoresistance 
measurements,  as  described  in  the  Supplementary  Section 
‘Ferromagnetic  contacts’.  The  Fe  contacts  exhibited  square 
magnetization  versus  applied  field  loops  with  a  coercive  field  of 
62  +  5  Oe.  The  Co  contacts  exhibited  slightly  skewed  loops  with 
much  higher  coercive  fields  of  150  +  5  Oe. 

To  probe  the  current-induced  spin  polarization  of  the  Tl  surface 
states,  measurements  were  performed  with  an  Fe/Al203  contact  as  a 
detector  on  a  45-nm-thick  Bi2Se3  film,  using  the  contact  layout 
shown  in  Fig.  2b.  A  fixed  bias  current  flows  between  two  other  con¬ 
tacts  on  either  side  of  and  co-linear  with  the  detector.  The  voltage  at 
the  detector  was  recorded  with  a  high-impedance  voltmeter 
(>1  GQ)  as  a  function  of  an  in-plane  magnetic  field  applied 
either  parallel  or  orthogonal  to  the  electron  bias  current  direction 
in  the  Tl.  No  current  flows  through  the  detector  contact.  A  linear 
background  was  subtracted  and  the  data  centred  on  the  vertical 
axis.  An  example  of  the  raw  data  and  background  subtraction  are 


shown  in  the  Supplementary  Section  ‘Raw  data  and  background 
subtraction’.  Electrons  flow  from  left  to  right  along  the  +x  axis 
for  a  negative  bias  current  and  a  positive  magnetic  field  points  in 
either  the  +x  or  +y  direction. 

When  the  magnetic  field  is  applied  orthogonal  to  the  current 
direction,  the  detector  contact  magnetization  is  switched  to  be 
either  parallel  or  antiparallel  to  the  current-induced  spin  in  the  Tl 
surface  (Fig.  lc)  and  the  detector  voltage  exhibits  clear  hysteretic  be¬ 
haviour,  as  shown  in  Fig.  3a  for  contact  Fe2.  The  relative  orientation 
of  Ie,  induced  Tl  spin  polarization  s  and  contact  magnetization  M 
are  shown  in  the  insets.  At  large  positive  fields,  the  magnetization 
is  parallel  to  the  Tl  surface-state  spin  polarization,  and  a  constant 
negative  voltage  is  measured  as  the  magnetic  field  decreases  to 
zero.  As  the  magnetic  field  changes  to  small  negative  values  (red 
trace),  an  abrupt  increase  in  detector  voltage  is  observed  at  the  coer¬ 
cive  field  (  —  65  Oe)  of  the  detector,  that  is,  when  the  detector  mag¬ 
netization  reverses  to  be  antiparallel  with  the  Tl  surface-state  spin. 
Further  ramping  the  magnetic  field  to  —800  Oe  produces  little 
change  in  the  detector  voltage,  and  the  overall  scan  exhibits  a 
single  step-like  behaviour.  When  the  field  sweep  direction  is 
reversed  so  that  it  increases  from  negative  to  positive  values  (black 
trace),  the  detector  voltage  is  constant  until  the  positive  coercive 
field  of  +65  Oe  is  reached.  At  this  point,  the  voltage  abruptly 
decreases  as  the  contact  magnetization  again  switches  and 
becomes  parallel  with  the  spin  orientation  of  the  Tl  surface  state. 
Similar  step-like  behaviour  is  observed  at  an  adjacent  detector 
contact  labelled  Fel,  as  seen  in  Fig.  3b. 
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Figure  4  |  Tl  spin  polarization  detected  as  a  voltage  with  Fe/Al203  and  Co/MgO/graphene  contacts.  The  contact  layout  of  Fig.  2c  is  used  for  two  Bi2Se3 
film  thicknesses,  a,  Optical  image  of  one  of  the  devices  measured.  The  bond  pads  near  the  bottom  of  the  figure  have  dimensions  of  100  x  100  pm. 
b,  Magnetic  field  dependence  of  the  voltage  measured  using  Fe/Al203  tunnel  barrier  detectors  on  an  11-nm-thick  Bi2Se3  film  with  the  magnetization  collinear 
with  the  induced  Tl  spin  for  a  bias  current  of  -0.5  mA.  c-f,  As  in  b,  but  with  Co/MgO/graphene  detectors  on  a  23-nm-thick  Bi2Se3  film  for  bias  currents  of 
+1  mA  (c),  -1  mA  (d),  +2  mA  (e)  and  —2  mA  (f).  g,h,  Field  dependence  of  the  voltage  measured  with  the  Co  detector  magnetization  orthogonal  to  the 
induced  Tl  spin  for  bias  currents  of  +2  mA  and  -2  mA,  respectively.  The  detector  contacts  have  areas  of  80  x  80  pm  and  all  data  were  obtained  at  T=  8  K. 
Black  traces  correspond  to  positive  magnetic  field  sweeps  (from  negative  to  positive  values)  and  red  traces  correspond  to  negative  field  sweeps.  The  y  axis 
label,  V(pV),  is  the  same  for  all  panels. 


When  the  direction  of  the  bias  current  is  reversed,  the  orientation 
of  the  Tl  spin  is  also  reversed  and  the  field  dependence  of  the 
FM  detector  voltage  flips  about  the  zero-voltage  axis,  as  seen 
in  Fig.  3c,d.  The  field  dependence  exhibits  the  same  step-like  behav¬ 
iour  as  observed  in  Fig.  3a,b,  and  for  this  current  direction 
closely  resembles  the  magnetic  hysteresis  loop  of  the  Fe  detector 
contact  itself.  The  detector  voltage  often  exhibits  a  local  peak  at 
the  coercive  field  as  the  contact  magnetization  reverses,  which 
we  tentatively  attribute  to  reversal  of  magnetic  domains  in  the 
detector  (which  is  not  single  domain).  The  voltage  on  a  non-mag- 
netic  detector  exhibits  no  step-like  behaviour  (see  Supplementary 
Section  ‘Control  experiments’). 

Transport  measurements  were  also  performed  with  the  applied 
field  and  detector  contact  magnetization  parallel/antiparallel  to 
the  current  direction,  and  these  data  are  summarized  in  Fig.  3e,f. 
In  this  configuration,  the  Tl  surface-state  spin  is  now  orthogonal 
to  the  detector  contact  magnetization,  and  a  completely  different 
field  dependence  is  observed.  No  abrupt  step-like  change  or  hyster¬ 
esis  in  the  detector  voltage  occurs  as  the  applied  field  is  swept 
through  the  contact’s  coercive  field  and  the  contact  magnetization 
is  reversed.  Instead,  a  smoothly  varying  voltage  is  observed, 
typical  of  a  pseudo-linear  magnetoresistance  obtained  with 
non-magnetic  contacts. 

A  second  set  of  devices  was  prepared  with  Fe/Al203  tunnel 
barrier  contacts  on  an  1 1  nm  Bi2Se3  epilayer  on  a  sapphire  substrate 
using  the  contact  layout  of  Fig.  2c.  To  improve  signal  quality,  a  non¬ 
magnetic  contact  immediately  adjacent  to  the  ferromagnetic  detec¬ 
tor  was  used  as  the  voltage  reference  to  directly  subtract  non-mag- 
netization-dependent  background  contributions;  this  results  in  a 
substantial  improvement  to  the  signal  quality  and  reduction  in  the 
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background  voltage.  An  optical  image  of  one  of  the  devices 
measured  is  shown  in  Fig.  4a,  and  representative  data  for  —0.5  mA 
bias  current  are  shown  in  Fig.  4b.  The  dependence  of  the  Fe  detector 
voltage  on  the  magnetic  field  and  direction  of  the  bias  current  is  the 
same  as  described  for  Fig.  3.  Note  that  even  though  a  smaller  unpo¬ 
larized  current  is  applied,  the  measured  spin  voltage  is  comparable 
to  that  for  the  45-nm-thick  Bi2Se3  film  of  Fig.  3. 

The  correlation  between  detector  magnetization  and  current- 
induced  spin  polarization  of  the  spin-momentum  locked  Tl 
surface  bands  is  more  easily  followed  using  a  ferromagnetic  detector 
metal  with  a  larger  coercive  field  and  a  skewed  magnetic  hysteresis 
loop,  so  that  reversal  of  the  contact  magnetization  occurs  over  a 
larger  range  of  applied  field.  To  this  end,  a  third  set  of  devices 
was  prepared  with  ferromagnetic  contacts  consisting  of  polycrystal¬ 
line  Co  films  on  2-nm-thickMgO  deposited  onto  a  graphene  mono- 
layer  previously  mechanically  transferred  onto  a  23-nm-thick  Bi2Se3 
film.  The  contact  geometry  of  Fig.  2c  was  used  and  the  graphene  was 
removed  between  the  contacts. 

These  data  are  summarized  in  Fig.  4c-h.  The  dependence  of  the 
Co  detector  voltage  on  the  magnetic  field  and  direction  of  the  Tl 
bias  current  are  as  described  previously,  demonstrating  that  the 
spin-voltage  behaviours  measured  are  reproduced  with  a  completely 
different  detector.  The  same  step-like  behaviour  in  the  measured 
spin  voltage  is  observed  as  the  contact  magnetization  switches  at 
Co  film  coercive  fields  of  +150  Oe  for  bias  currents  of  +1  and 
+2  mA  (Fig.  4c, e).  When  the  Tl  bias  current  is  reversed,  the  spin 
voltage  changes  sign  (Fig.  4d,f).  The  larger  coercive  field  of  the 
Co  and  the  slower  reversal  and  saturation  of  the  magnetization 
are  reflected  in  the  width  of  the  hysteretic  loop  observed  in  the 
detector  voltage.  This  is  especially  evident  when  the  magnetic 
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Figure  5  |  Dependence  of  the  ferromagnetic  detector  voltage  on  bias  current,  thickness  and  temperature,  a,  Height  of  the  voltage  hysteresis  curve 
AV  =  VCM^—Vi—M)  for  applied  fields  greater  than  the  contact  coercive  field  as  a  function  of  bias  current  for  samples  with  80  x  80  pm  Co/MgO/graphene 
detector  contacts,  b,  AR  =  \AV/Ib\  as  a  function  of  the  Bi2Se3  film  thickness,  where  lb  is  the  bias  current  applied  through  the  Bi2Se3.  Dashed  line  shows  an 
inverse  thickness  dependence  scaled  vertically  to  fit  the  data.  Each  data  point  corresponds  to  a  different  bias  current  for  a  device  on  each  of  the  three  Bi2Se3 
film  thicknesses  used  here.  Green  diamond  and  blue  triangle  symbols  are  for  Bi2Se3  films  with  Fe/Al203  contacts,  while  the  red  circles  indicate  Bi2Se3  films 
with  Co/MgO/graphene  contacts,  c,  Height  of  the  voltage  hysteresis  curve  |AV|  as  a  function  of  temperature  for  samples  with  Co/MgO/graphene  detector 
contacts  on  a  23  nm  Bi2Se3  film  and  2  mA  bias  current.  The  detector  voltage  is  measured  with  the  magnetization  col  linear  with  the  induced  Tl  spin, 
d,  Detector  voltage  versus  magnetic  field  at  125  K  for  the  Co/MgO/graphene/23  nm  Bi2Se3  film  and  +2  mA  bias  current.  The  black  trace  corresponds  to  a 
positive  magnetic  field  sweep  (from  negative  to  positive  values)  and  the  red  trace  corresponds  to  a  negative  field  sweep.  The  error  bars  in  panels  a  and  c 
represent  the  uncertainty  in  determining  the  height  of  the  voltage  hysteresis  curve  arising  from  the  signal-to-noise  ratio  and  the  background. 


field  is  applied  along  the  bias  current  direction,  as  shown  in 
Fig.  4g,h.  At  high  applied  field,  the  magnetization  is  orthogonal  to 
the  current-induced  spin  polarization  and  the  projection  of  the 
surface-state  spin  onto  the  detector  magnetization  is  zero.  As  the 
applied  field  is  reversed,  the  magnetization  rotates  through  a  parallel 
(or  antiparallel)  orientation  with  the  spin,  producing  a  small  peak  in 
the  detector  voltage,  and  then  again  becomes  orthogonal  to  the 
induced  spin.  The  linear  background  evident  in  Fig.  3e,f,  attributed 
to  conventional  magnetoresistance,  is  absent  here  because  the  non¬ 
magnetic  contact  adjacent  to  the  detector  is  used  as  reference  for  the 
voltage  measurement. 

The  change  in  detector  voltage,  Ay  =  V(M)  —  V(—M),  measured 
as  the  magnetization  reverses,  increases  linearly  with  bias  current  Ih 
flowing  in  the  Bi2Se3  film,  as  shown  in  Fig.  5a  for  the  Co-based 
detector  data  of  Fig.  4.  The  value  |AV//b  |  =  A R  scales  inversely 
with  the  thickness  of  the  Bi2Se3  film  (Fig.  5b),  demonstrating  that 
the  detector  voltage  measured  derives  from  the  TI  surface  states 
and  not  the  bulk.  The  bias  current  applied  between  the  non-mag- 
netic  ohmic  contacts  is  shunted  by  the  bulk  and  only  that  fraction 
flowing  in  the  surface  layer  contributes  to  spin  momentum 
locking  in  the  TI  surface  state  to  generate  a  spin  voltage.  If  the 


sample  consists  of  N  quintuple  layers  (QL,  the  basic  five-atom- 
layer  unit  Se-Bi-Se-Bi-Se,  ~1  nm  thick  in  this  naturally  layered 
compound)  and  each  QL  has  equal  conductivity,  then  the  bias 
current  through  the  surface  layer  is  1/M  and  the  spin  voltage 
should  scale  accordingly.  This  behaviour  is  clearly  observed 
in  Fig.  5b. 

These  observations  for  the  two  ferromagnetic  contact  structures 
underscore  the  reproducibility  of  the  data  and  the  fact  that  these 
behaviours  are  due  to  bias  current-induced  spin  in  the  TI  Dirac 
surface  states  and  are  independent  of  the  details  of  the  contact 
itself.  The  flow  of  current  in  the  TI  surface  band  induces  a  net 
spin  polarization  aligned  in-plane  and  orthogonal  to  the  current. 
The  strong  correlation  of  the  measured  ferromagnetic  detector 
voltage  with  both  the  amplitude  and  the  relative  orientation  of  the 
bias  current-induced  spin  to  the  detector  magnetization  provide 
clear  and  compelling  evidence  for  spin-dependent  detection. 

This  also  excludes  possible  contributions  from  local  Hall  effects 
resulting  from  perpendicular  components  of  the  fringe  field  at  the 
edges  or  domain  boundaries  of  the  ferromagnetic  detector 
contact,  or  other  Lorentz  force-based  effects,  because  any  such 
fringe  fields  are  present  regardless  of  the  orientation  of  the 
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in-plane  magnetization  and  direction  of  current  flow.  The  voltages 
on  non-magnetic  detector  contacts  and  magnetic  ohmic  (rather 
than  tunnel  barrier)  detector  contacts  show  no  such  field  dependent 
switching  (see  Supplementary  Section  ‘Control  experiments’), 
excluding  contributions  from  other  mechanisms.  We  also  discuss 
and  exclude  contributions  from  anisotropic  magnetoresistance 
and  anomalous  Hall  effects,  which  might  arise  from  current  shunt¬ 
ing  by  the  detector  contact  (see  Supplementary  Section  ‘Anisotropic 
magnetoresistance  and  anomalous  Hall  effects;  and  Additional 
control  experiments’).  Thus,  the  voltage  on  our  ferromagnetic 
tunnel  contact  directly  reflects  the  existence  and  orientation  of  a 
current-induced  spin  polarization  of  the  TI  Dirac  surface  bands. 

The  hysteretic  behaviour  of  the  detector  voltage  decreases  with 
increasing  temperature  and  disappears  by  150  K.  The  height  of 
the  hysteresis  curve  | V(M)  —  V(— M)|  for  applied  fields  greater 
than  the  contact  coercive  field  is  plotted  as  a  function  of  temperature 
in  Fig.  5c.  The  data  at  125  K  are  shown  in  Fig.  5d.  Spin-resolved 
photoemission  has  mapped  the  Bi2Se3  surface  states  at  temperatures 
as  high  as  300  K  (refs  14,23).  This  discrepancy  may  be  due  to  effects 
induced  by  the  ambient  measurement  environment  used  for  our 
work  rather  than  the  ultrahigh-vacuum  environment  required  for 
photoemission,  or  to  mechanisms  that  are  transport  specific. 
Because  our  samples  are  strongly  n-type  (~1  x  1019  cm  3),  we  ten¬ 
tatively  attribute  the  temperature  limit  of  150  K  to  thermal  contri¬ 
butions  from  the  bulk  that  dilute  and  overwhelm  the  surface  band 
contribution.  A  more  meaningful  assessment  of  temperature  depen¬ 
dence  will  require  better  control  of  the  carrier  density  either  by 
electrostatic  gating  or  improvements  in  sample  fabrication. 

The  use  of  a  ferromagnetic  contact  in  a  multi-terminal  configur¬ 
ation  to  detect  a  net  spin  polarization  produced  in  an  adjacent 
non-magnetic  medium  by  some  form  of  spin  injection  is  well  estab¬ 
lished28  and  has  been  applied  to  a  variety  of  materials  including 
metals,  semiconductors  and  carbon  nanostructures30,31.  The  ferro¬ 
magnetic  contact  is  typically  used  as  a  voltage  probe  that  draws 
no  current  (the  voltage  is  proportional  to  the  projection  of  the 
spin  polarization  onto  the  detector  magnetization).  Such  an 
approach  has  been  quantitatively  treated  by  Hong  and  colleagues 
for  the  specific  case  of  current-induced  spin  polarization  in  a  TI 
in  both  diffusive  and  ballistic  regimes29.  They  show  that  the  ferro¬ 
magnetic  detector  voltages  V(M)  we  measure  experimentally  (Figs  3 
and  4)  are  directly  related  to  the  current-induced  spin  polarization 
by  [V(M)  — V(  — M)]  =  JbRBPFM(p-Mu),  where  Ib  is  the  (hole)  bias 
current  in  the  +x  direction,  RB  is  the  ballistic  resistance  of  the 
channel  and  Am  is  the  transport  spin  polarization  of  the  FM  detec¬ 
tor  metal.  Mu  is  a  unit  vector  along  the  detector  magnetization  M, 
and  p  is  the  degree  of  spin  polarization  induced  per  unit  bias 
current  by  both  spin-momentum  locking  in  TI  Dirac  surface 
states17,18  as  well  as  Rashba  spin-orbit  coupling.  The  Dirac  surface 
state  and  Rashba  contributions  are  opposite  in  sign. 

Our  measured  values  [V(M)  —  V(  — M)]  exhibit  the  essential 
characteristics  expected  from  this  equation:  they  scale  linearly 
with  bias  current  and  exhibit  the  angular  dependence  given  by 
p-Mu.  The  sign  of  our  measured  values  [V(M)  —  V{— M)]  indicates 
that  the  Dirac  surface  state  contribution  dominates,  as  expected 
from  theory18  and  explicit  model  calculations29.  We  discuss  the 
sign  of  the  voltage  expected  for  spin-momentum  locking  and 
Rashba  in  the  Supplementary  Section  ‘Two  dimensional  electron 
gas  states  and  potential  Rashba  contribution’  and  show  that  the 
Rashba-induced  polarization  is  expected  to  be  small,  even  for 
strongly  enhanced  spin-orbit  coupling.  The  temperature  depen¬ 
dence  we  observe  in  Fig.  5c  is  also  significantly  stronger  than  that 
reported  for  Rashba-induced  polarization32.  A  quantitative  exper¬ 
imental  determination  of  p  is  more  elusive,  because  the  bias 
current  has  parallel  bulk  and  surface-state  conduction  paths  and 
only  that  fraction  flowing  in  the  surface  channel  contributes  to  the 
spin  polarization  arising  from  Dirac  states.  We  estimate  |p|  «  0.2 
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from  the  data  of  Fig.  4b  or  e  based  on  the  following  assumptions: 
PFM(Fe,  Co)  «  0.4,  kF«0.15A  ',  only  bias  current  flowing  in  a 
1-nm-thick  surface  channel  contributes  to  the  spin  polarization, 
and  the  conductivities  of  the  surface  and  bulk  conduction  channels 
are  equal.  This  value  is  reasonably  consistent  with  theory18,29,  but 
subject  to  the  uncertainties  in  bias  current  distribution  already 
noted.  It  represents  a  lower  bound  for  the  polarization  produced 
by  spin-momentum  locking,  because  in  our  highly  n-doped 
samples,  £F  lies  well  above  the  conduction  band-edge  and  contri¬ 
butions  from  Rashba  spin-split  states  with  opposite  net  spin  polar¬ 
ization  will  reduce  the  voltage  we  measure  and  the  corresponding 
estimate  of  TI  surface-state  polarization  (fig.  4  of  ref.  29). 

In  conclusion,  the  intrinsic  sensitivity  to  spin  polarization 
afforded  by  a  ferromagnetic  contact  enables  one  to  distinguish  the 
current-induced  TI  surface  band  spin  polarization  from  an  unpolar¬ 
ized  background.  The  magnitude  and  temperature  dependence  of 
the  spin  voltage  we  observe  is  probably  compromised  by  bulk 
(non-TI  surface)  states,  which  provide  a  parallel  conduction  path 
and  a  large  unpolarized  electron  background.  We  anticipate  that 
improvements  in  Bi2Se3  growth  will  lead  to  significantly  lower  unin¬ 
tentional  doping  levels,  reducing  parasitic  conduction  pathways  and 
enhancing  the  effects  we  observe.  Gating  such  samples  to  position 
£f  between  the  Dirac  point  and  the  bulk  conduction  band  is 
expected  to  significantly  increase  the  spin  voltage  and  measured 
spin  polarization  and  will  enable  further  insight  into  the  behaviour 
of  the  TI  surface  states  and  spin-momentum  locking. 

Methods 

The  growth  of  Bi2Se3  films  was  carried  out  on  epitaxial  graphene/SiC  (0001)  and 
Al203(0001)  substrates  in  an  ultrahigh-vacuum  (UHV)  system  (base  pressure 
~  1  x  10- 10  torr)  that  integrates  two  MBE  chambers  and  a  low-temperature 
(5-300  K)  scanning  tunnelling  microscope  (STM)33,34.  As-received  nitrogen- doped 
6H-SiC(0001)  substrates  were  first  etched  in  a  H2/Ar  atmosphere  in  a  separate 
chamber  at  <1,500  °C  to  remove  polishing  damage.  They  were  then  transferred  to 
the  MBE  system  and  annealed  in  UHV  in  a  Si  flux  (<0.1  monolayers  min-1)  at 
950  °C  to  produce  a  (3  x  3)  reconstructed  surface  and  further  annealed  at 
temperatures  of  1000-1300  °C  without  Si  flux  to  grow  epitaxial  graphene35.  For  the 
growth  of  Bi2Se3,  the  substrate  was  held  at  275-325  °C  and  Bi  and  Se  were  supplied 
via  separate  Knudsen  cells  at  460  and  250  °C,  respectively.  For  Bi2Se3  growth  on 
Al203(0001)  substrates,  a  two-step  process  was  used36:  2-3  QL  of  Bi2Se3  were 
first  deposited  at  a  reduced  temperature  of  100  °C  and  the  substrate  temperature 
was  then  slowly  raised  to  300  °C  where  the  rest  of  the  film  was  deposited.  In  situ 
STM  imaging  was  used  to  monitor  the  surface  morphology  and  electronic 
structure  and  to  ensure  optimal  layer-by-layer  spiral  growth.  The  as-grown  film 
exhibited  a  Dirac  point  250  meV  below  the  Fermi  level  in  tunnelling  spectra, 
indicating  the  n-type  doping  typically  seen  in  MBE-grown  materials  due  to  Se 
vacancies.  The  character  of  the  topologically  protected  surface  states  has  been  shown 
to  be  insensitive  to  morphological  defects  such  as  grain  boundaries34. 

Two  types  of  ferromagnetic  detector  contacts  were  used  as  spin-sensitive  probes. 
Such  contacts  have  successfully  enabled  electrical  detection  of  spin  current  and 
accumulation  in  both  semiconductors  and  metals28,30,37-42. 

Fe/Al203  contacts  were  formed  on  the  air-exposed  surface  of  11 -nm-  and 
45-nm-thick  Bi2Se3  films  in  a  separate  MBE  system  as  follows39.  A  0.7  nm  layer  of 
polycrystalline  A1  was  first  deposited  by  MBE  and  then  oxidized  in  200  torr  02  for 
20  min  in  the  presence  of  ultraviolet  light  in  the  load-lock  chamber.  This  step  was 
then  repeated  for  a  total  A1203  thickness  of  2  nm.  The  sample  was  then  transferred 
under  UHV  to  an  interconnected  metals  MBE  chamber,  where  20  nm  of 
polycrystalline  Fe  was  deposited  at  room  temperature  from  a  Knudsen  cell.  An  oxide 
tunnel  barrier  was  used  to  prevent  potential  direct  magnetic  exchange  with  the  TI 
surface  states,  interdiffusion  (which  may  result  in  undesirable  alloy  formation 
and  spin  scattering)  and  to  better  sample  the  TI  surface  state  at  the  interface 
rather  than  the  bulk.  The  samples  were  processed  into  the  device  structures 
illustrated  in  Fig.  2  to  enable  transport  measurements.  Standard  photolithography 
and  chemical  etching  methods  were  used  to  define  the  Fe  contacts,  which  ranged  in 
size  from  10  x  10  pm2  to  80  x  80  |xm2,  with  adjacent  contact  separation  ranging 
from  45  to  200  |xm.  Ion  milling  was  used  to  pattern  the  Bi2Se3  mesa.  Large 
Ti/Au  contacts  were  deposited  by  liftoff  in  an  electron-beam  evaporator  as  non¬ 
magnetic  reference  contacts  and  bias  current  leads.  The  Fe  contacts  were  capped 
with  10  nm  Ti/100  nm  Au,  and  bond  pads  for  wire-bonded  electrical  connections 
were  electrically  isolated  from  the  substrate  using  100  nm  of  Si3N4. 

Another  set  of  samples  were  prepared  with  Co/MgO/graphene  as  the 
ferromagnetic  detector  contacts  on  a  23-nm-thick  Bi2Se3  film.  We  have  previously 
shown  that  single-layer  graphene  forms  an  effective  tunnel  barrier  on  both  metal41 
and  semiconductor42  surfaces,  which  prevents  oxidation  and  interdiffusion43.  A 
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single  graphene  layer  grown  by  chemical  vapour  deposition44  was  transferred  onto 
the  air- exposed  Bi2Se3  surface  and  then  patterned  so  that  it  was  only  located  under 
the  detector  contacts  (both  magnetic  and  non-magnetic).  A  2  nm  MgO  film  was 
then  deposited  by  electron-beam  evaporation  following  the  deposition  of  a  0.1  nm 
Ti  seed  layer.  A  15  nm  Co  layer  was  immediately  deposited  in  the  same 
evaporation  system  without  breaking  vacuum.  Similar  photolithography  steps  as 
described  above  were  used  to  produce  the  device  structures  illustrated  in  Fig.  2,  with 
adjacent  rows  of  magnetic  and  non-magnetic  (reference)  detector  contacts.  Separate 
magnetometry  measurements  confirmed  that  the  ferromagnetic  contact 
magnetization  lay  in-plane  with  square  hysteresis  loops  for  the  Fe  and  slightly 
rounded  loops  for  the  Co,  as  expected  for  these  15-  to  20-nm-thick  films  and  their 
respective  magnetocrystalline  anisotropies.  The  contact  resistances  of  the  Au/Ti 
reference  and  magnetic  detector  contacts  were  ~  10  IT  and  103-104  fl,  respectively. 

Transport  measurements  were  performed  in  a  closed-cycle  cryostat  equipped 
with  an  electromagnet  (4-300  K,  +1,000  Oe).  A  bias  current  of  0.5  to  3  mA  was 
applied  through  a  pair  of  contacts  on  either  side  of  and  collinear  with  a  selected 
detector  contact,  and  the  voltage  on  the  detector  contact  was  recorded  as  a  function 
of  the  in-plane  magnetic  field  applied  either  parallel  or  orthogonal  to  the  electron 
bias  current  direction  in  the  TI. 
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